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АВЗТКАСТ 
The stresses near the root of an F8U-3 airplane wing loaded in 

torsion were found using the matrix-force method of analysis. The wing 
was represented by a simplified, idealized model consisting of axial- 
load carrying flange bars and shear carrying cover plates. The theoret- 
ical stresses in the idealized model were found by using the IBM 360/67 
computer to solve for the redundant loads. These stresses show favor- 
able comparison with the experimental stresses determined from strain 
gages mounted on the cover skin of the wing. This study shows that the 


matrix-force method is applicable to thick-skinned wings. 
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INTRODUCTION 


The purpose of this study is to predict the stresses im an F8U-3 
airplane wing loaded im torsion using the matrix-force method of analysis. 
The stresses determined analytically are compared with the stresses ob- 
tained experimentally in order to determine the accuracy and consistency 
of the results. 

The primary structure of the F8U-3 wing center section consists of 
a simple multicell torque box, as shown in Figure 1. Spanwise bending 
is carried by the relatively thick skins which are fabricated from 
7079-T6 aluminum alloy plate. ` The skins are machined to varying thick- 
nesses both spanwise and chordwise in order to obtain an optimum distri- 
bution of material. The skins are continuous from the wing fold rib 
inboard to the wing centerline and from the front beam to the rear beam. 
It is this area, between wing fold ribs, port and starboard, and Front 
and Rear Beams, fore and aft, which defines the wing center section, 
the area of interest in this study. The sweep of the wing is 42° at 
the quarter chord. The wing thickness varies from 4% at the tip to 
5.01% at the ποθ. It is this combination of a high degree of sweep 
and low thickness ratio which raises problems with structural analysis 
of the wing. Due to bending-torsion interaction, or sweep-coupling, 
structural analysis by engineering beam theory methods runs into diffi- 
culty when analyzing stresses near the wing root. 

This thesis describes the matrix-force method of analysis for 
determining the stress distribution in such a wing. Among the advan- 
tages of this method are: 

(1) It results in a consistent set of internal stresses and external 


deflections. 
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FIG. ! 
εθυ-2 WING CENTER SECTION 


(2) The structure can be altered and re-analyzed with only minor 

changes in the computer data cards. 

(3) The "Bookkeeping" can be simplified and organized in a logical 

manner. 

The matrix-force method of analysis uses forces as the unknowns and 
treats the structure as an assembly of elastic components. The method, 
as described in Reference 1, is divided into three parts: 

A. Idealize the structure 
1. Replace the actual structure with a mathematical model. 


2. Make the model statically determinate by cutting redundant 
members. 


3. Write the equilibrium equations that determine the load distribu- 
tion in the structure 


B. Calculate and tabulate the flexibilities of the individual members 
comprising the structure. 
C. Perform matrix operations. 

In this report, the wing is "cut" outboard of the intermediate 
rib Y, = 81.98), and it is assumed that the stress distribution here 
can be accurately calculated using the engineering theory of torsion, 
7 = = 24. The dihedral of the wing is neglected, and the wing is 
assumed to be supported on four rigid supports located at A, B, C, 
and D in Figure 2. The aft supports are pins which can support fore 
and aft and vertical loads. The forward supports can support vertical 
loads only. The only loading condition considered in this study is 
symmetrical wing torsion. From symmetry, the reactions at all four 
support points are determined by equilibrium as functions of the applied 


loads. 
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This study was conducted at the Naval Postraduate School, Monterey, 
California, during 1967. The author is indebted to Professor C. H. Kahr 
who conceived the need for the investigation and who rendered valuable 
assistance and encouragement throughout the project. The author expresses 
gratitude to Aeronautical Laboratory Technicians R. A. Besel and T. B. 
Dunton whose assistance and cooperation were essential to the completion 


of the study. 
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The model selected to represent the wing structure consisted of 
flat plates which carry only shear and cylindrical bars which carry 
only axial load. This model was selected for compatibility with the 
previous work done in Reference 2. Analysis of the model can be as 
accurate as the matrix-force method and arithmetic permit; however, 
the correlation between theoretical and experimental stresses will 
depend upon how closely the model represents the actual structure. 

This study is not an attempt to define the best model for the analysis, 
but to calculate stresses in what was tacitly assumed to be an adequate 
model. Most of the examples of similar studies currently available in 
the literature deal with thin-skinned shells. There is little evidence 
on whether or not such models are suitable for the thick-skinned wing 
plates used in this study. Ап important factor in constructing the 
idealized model is the change of the wing geometry to model geometry. 
The model consists of bars and flat shear carrying plates. It is ex- 
tremely important to insure that these plates are, in fact, mathemat- 
ically flat. For example, the depths of the ribs at WS 25.095 and WS 
81.98 were such that the four points AFGH in Figure 1 did not lie in a 
plane. Therefore,panel AFGH is warped and has panel edge shears which 
have components perpendicular to plane AFGH. These shears are not 
accounted for in the equilibrium equations written for a flat panel 
and will produce errors in the solution. Reference 3 shows that for 
a similar wing such discrepancies yielded determinate member loads 
which were out of equilibrium as much as 5 to 10%. After altering 
the beam depths to make the panels flat, this trouble disappeared. The 


alteration required a change of less than 1% of each dimension. 


After determining model geometry from the adjusted dimensions of 
the wing, the cross sectional area of the bars and the thickness of the 
shear panels were calculated. Each bar area of the model consists of 
the area of the beam flange plus one half the skin area between adjacent 
beams. The beam web areas were not included, using the methods of 
Reference 4. The effect of combining the axial load carrying skin with 
the beam flange to form the model's cylindrical bar is to remove Poisson's 
effect from the analysis. This means that any chordwise contraction of 
the cover skins is neglected. The wing skin thickness varies in the 
chordwise direction in section ABCD and in both spanwise and chordwise 
directions in section ABEF. An average skin thickness between beams 
was determined and used as the shear plate thickness. Using constant 
thickness shear plates in place of tapered skin simplifies the analysis 
but introduces possible sources of error. Model geometry and areas are 
listed in Tables I and II. Diagrams of cross sectional geometry are shown 
in Figure 3, 4, and 5. An exploded view of the outer section of the 
idealized model from ye - 25.095 to Үс = 81.98 is shown in Figure 6. 
An exploded view of the intermediate rib at Y, = 81.98 is shown in 
Figure 7. 

2. Make the model statically determinate by cutting redundant members. 

The most complex beam that is truly determinate, so that the flange 
and shear loads can be determined from equilibrium equations, is the 


3-flange, 3 web beam shown in Figure 8. 


Fig 8. 


DETERMINATE BEAM 
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In this beam, loads can be obtained from equilibrium only, without con- 
sidering flange areas, material stiffness, or deformation of the structure. 
In order to make the idealized model of this study a determinate structure, 
it was mecessary to cut the 5 interior webs and 11 of the flanges, re- 


sulting in the structure shown in Figure 9. 
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Fig. 9 


CROSS SECTION OF DETERMINATE SYSTEM 


The number of redundancies arising from the geometry of a box structure 
has been given by Reference 4 and others as B + N - 4 where B is the 
number of longitudinal effective flanges continuous across each junction 
and N is the number of closed cells stiffened by ribs at tip end of bay. 
For the outer section of this model, B = 14 and N = 6,which results in 
16 degrees of redundancy. This agrees with the 16 members which were 
cut to provide the determinate structure of Figure 7. 

3. Write the equilibrium equations that determine the load distri- 
bution in the structure. 

Considering the outer wing section from station 25.095 to 81.98, 
there are the following 49 unknown member loads. 


12 skin panels 
14 beam flanges 


7 beam webs 

10 rib flange loads at station 81.98 
6 rib webs at station 81.98 

49 total unknowns 


The redundant forces come from the following members. 


1 F.B. , upper flange 
2 FR.I.B., both flanges 
2 F.I.B. , both flanges 
2 C.I.B. , both flanges 
2 A.I.B. , both flanges 
2 R.1.B. , both flanges 
5 webs 8 - 12 


redundant forces 


ON 


Thus, there are &9 - 16 or 33 equations of equilibrium required to 
characterize the model. These equations are: 
a) Six equations of equilibrium of free body of bay between stations 


Y = 25.095 and 81.98 


w 
(1) P =O 
2) P,=0 
E 
(4) M = 0 
6) M. -0 
(6) M = 0 


b) Eleven equations of equilibrium of the redundant beam flanges, elements 
16, 18, 29; 22, 24, 206 3485 39$ 35, 37, and 39. 

c) Ten equations of forces and moments, F. - 0, ы” = 0 of the 
intermediate beam segments of the rib at Y, = 81.98 


d) Six equations of equilibrium of the rib posts on the beam ends at 


Y = 81.98, elements 53 through 58. 


The addition of the center section Y. O to 25.095 can be thought 
of as an extention of the horizontal flanges. The forces in the X and Z 


directions are reacted by the horizontal and vertical supports at Y, = 


25.095, but the Y forces carry through to the centerline, Y, = 0. 
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Although the center section is 16° redundant (the same as the там 
board section), these are not new redundancies but merely components 
of the forces in the outer section. Further investigation shows that for 
symmetrical wing loading, the bending meoment, М, is constant from 
id. - 25.095 to Y = 25.095 and that torsional and vertical shear are 
zero. The axial load on each flange at E - 25.095 is the same as the 
load on each flange at Y. = 25.095. Hence there is no shear on any 
panel or web in the center section. (It was later shown that due to 
a slightly unsymmetrical load, there are, in fact, shearing stresses 
in the center section panels.) Since the addition of the center section 
results in no new redundancies and no new external loads, there can be 
no more than the 33 independent equilibrium equations already proposed 
for the outer section. Since there are more than 33 structural elements 
involved, there are several possible combinations of 33 independent 
equations. Those chosen were considered to be those most likely to 
have the larger stress concentrations. 
B. Calculate and tabulate the flexibilities of the individual members 
comprising the structure. 

As described in Section C, the internal strain energy of the wing 
may be written as: 

a, 
AL ml dan. Е Ие, = (1) 
а а 

Тһе elements in the F matrix define the elasticity of the structure. 
If skin thicknesses or flange areas are changed, only the elements of 
the F matrix will be affected. If the geometry of the model is changed, 
the equations of equilibrium will also be affected. Multiplying out 


equation 1, 
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ZEU = ТЕ + a + За, ғара; "2858.05. τ. (2) 
Each term of Equation 2 can be regarded as a small increment of strain 
energy AU associated with the load a, ога, . а Since the load a, 
may be applied to more than one structural element, the corresponding 
ΠΗ term may be made up of contributions from several structural elements. 
The F matrix must be a square symmetric matrix where в = της 
1. Flange Bars. 

The calculation of the F matrix elements for the flanges was done 


following the example of Reference 4. The energy in a uniform bar under 


a linearly varying axial force 15 


2 2 
/ 2, - Q, 
U = siz] Кэт" d »" (3) 


where 9, апа еке end forces as shown in Figure 10. 


t Ze 


TM 


Fig. 10 


FLANGE BAR 


Rewriting Equation 2 for this single element 
e 2 
РЕЙ = ФТ а "4; "m e 


therefore 


Q^ 
Ni < 


Aez + ах (5) 





жыз ο 
r на 2; А 





ο 


"s. 
Pn ή (6) 
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and fi, = 2, а Езу ο a 
ο 
(7) 


In the case of linearly tapered members, the flexibilities may be 
obtained as functions of the end area ratios. Reference 1 shows the 
flexibility influence coefficients modified by a function, Q, of the 
end areas which was obtained from a graph of Q versus А, / А the end 


For most of the flange elements, the generalized end loads, Ч; 


апа 4. were linear functions of the unkown member loads, or: 


9; = =, 9, r. Q nz О Ch) (9) 
9; = ба, d gm m бла, (10) 
where &. = the coefficient of the i'th member load applied to one end 


( 


of the bar, 

and Ё; = the coefficient of the i'th member load applied to the other 
end of the bar. 

Now taking each force acting on the ends of the flange individually, 


that is: 


7, 2 сс, а, (11) 
Ze = б, а, (12) 
and substituting into Equation 4, we find that 


2 


Z = a, (e, A, + 20,6, te + 6, on ) (13) 
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and from Equation 1, for a single force, а, 


2 
δα = a, Бр (14) 


Comparing Equations 13 and 14, 


2 2 
Б, шо, ⁄ + 2o¿ 6, + (8, f, (15) 


Similarly, for the other forces acting on the same flange, 


2 2. 
б. = ού, £, + ыз 2 В. 2, ж ©, Lia (16) 


and for the cross product term, 


БА ri od, ed a, f, τ В.В. ^; + #, (6, 4; + Bz, ) (17) 


The forces on element 14 are shown in Figure 11. 
For example, the flexibility influence coefficients for element 14 were 


determined as follows: 


ЕТЕ ear 
А, / 8 абаны — — 


Fig. 11 
FORCE DIAGRAM OF FLANGE BAR 


6 
ғ, 7 34, д, = 20.333 = #, 


LL = 
“Z - c аш 


,: 2728425 το. -/0.2530, 


Te 


. 68 Aa °@,, 


: 2 
т. , 220.233 (-70.953) * 2 x vo." 7 (C 70. 253)(0) +20.333(0) = 2439. 483 


, 


2 
2 20 7 20.333(- 3.23/) * 20.333(0Y(. 568) » 10.167 (- 10.953X. SEB) = 656.335 


29 


For element 14, there are 4 more coefficients to be found; Foo. 20? 
2 


21, 21” πο. 21? and Foo, 21: Тһе Р), 2 coefficient for element 14 18 


added to other F, 2 coefficients for elements 7, 39, 101, and 126 to 


make up the total F 


F 


2. 2 flexibility influence coefficent found in the F 
matrix. 
2. Symmetrical Shear Panels 


Reference 3 gives the energy of the symmetrical trapezoidal 


panel in terms of the shear force on the shorter parallel side, as: 








Ly (w, + W.) / м, — Wa = 
кр HA а? Іа ΤΟ. (18) 
= + 
¿EU | w, t 6 ( La ) Š (/+p) a, 
so that 
а 
= a U 
ae = Э 2 (19) 
2, 
may be written for Р. 1 of element 1 
Fis 14.335 (12. 460) / ( 
A? PESA ora) + t: 
46.572 (. 200) |6 
The forces on element 1 are shown in Figure 12. 
ЕС 
= 2 
d (030, 
Fig? MW 
FORCE DIAGRAM OF BEAM WEB 
Other ΕΙ 1 coefficients appear from elements 40 and 51 and add to pro- 
duce the final F found in the F Matrix. 


ls #1 
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Тһе shear panel edge loads were determined by assigning a unit 
shear force, a.» to the outboard edge of all panels and the forward 
edge of the rib webs. From the equilibrium of the free body, forces 


а» а, and a, may be found as functions of the unit force а.. 


d 
ЗБ 7 бі І / Wa 
ас аа. (20) 


2g" 33 14 / М. 

The geometry of the panels and the load ratios are shown іп Table III. 
3. Non-symmetrical shear panels. 

a) The cover skin panels have roughly the shape of a parallelogram, 

with the inclined sides skewed slightly. The maximum difference in 

inclination, 4-6 = 3° 34', occurs at the forward panels, elements 38 

and 15. These panels were treated in accordance with Reference 6 


which shows that the energy of the swept panel is given as 





ми [omar ga © 
where “ў з cof 6 + соё $ соё В + соб d (22) 


The forces on element 15 are shown in Figure 13. 








Fig. 13 
FORCE DIAGRAM OF COVER PANEL 

so that b 

F c дс 

n,n = 
да, 
may be written for Foo, 20 of element 15. , 
$2 / 32 (22/99) 516 (7» ^ 268 ) (3. 590) = 90.525 
| 7 
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` 
b) The cover skin panels introduce additional forces due to interaction 
between torsion and bending. Twist of a swept wing produces bending and 
vice-versa. This interaction was introduced by considering the stress 
distribution in the rhomboidal cover skin panels, as shown in Figure 12. 
These panels were assumed to be parallelograms in order to simplify the 
analysis. Consider the plate shown, loaded with a uniform shear force, 


a , and axial load, і... 


СТ u θα. d = =S 
— — 4 —— AR. 
Ёны ή == 7 
4 
£ 


= 2a tan A 





Fig. 14 
SWEEP COUPLING OF COVER PANEL 
The axial force, f. was assumed concentrated at the flanges which were 
made up mostly of skin material, so that the cover panels carry omis 
shear. However, it is seen in Figure 12 that the shear forces along 


the sides of the panel produce additional tensile stresses, f,', which 


t 
do work on the displacements produced by the primary axial stress, Е... 


For this panel, the strain energy can be written. 


2 Е 
у (23) 
CEU = 27, = + EL 


The third term in Equation 22 is due to the interaction between the 
axial and shearing forces and does not appear in an equation. written 

for a rectangular panel. This term causes the stress distribution in 
this model to show a root stress concentration. Reference 3 states that 
"Studies made on simple one cell swept box models with streamwise ribs, 
utilizing parallelogram shaped cover sheets, show that unless this cross- 


product term is introduced, there will be nothing in the analysis to 
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differentiate between a straight and swept wing, and there will be no 
i 3 ; 
root stress concentration factor'". For the purpose of calculating the 


energies in the flanges adjacent to each panel, this axial load, Е, 
3 


was split equally and added to the end of each flange. These coupling 
loads are shown as dashed vectors in Figure 6. 
C. Perform matrix operations. 

The matrix-force method, as described in References 1 and 3, 
expresses the energy of the entire structure as the successive product 


of three matrices. 
геа з (a) (A {a} (24) 


where a, ..... а аге the internal loads acting on the members. The 

F matrix is symmetrical and represents the elasticity of the structure. 
The internal loads were found by writing the equilibrium equations for 
the determinate structure and for each redundant member. The number 

of equations must equal the total number of unknown internal loads less 
the number of redundants. The equations are written so that all the 
independent external loads and redundants are collected on one side and 
non-redundant loads on the other side of the equality; 


a] fa,,) = [8] (e 12, (25) 


where [А] and [5] are matrices of the coefficients of the loads. 
Solving for the non-redundant loads, а >» 25 functions of the applied 


loads, P , and the redundants, a, » ме get: 


=l 


(a. - [le] eta.) - Mo 1a] as 


B9 


Then, by adding the trivial relation between the redundants 


(2,} = (г, (27) 
and expanding М to include Equation & results in 
(аў 7 [5] 2 12, (28) 
which may be transposed, yielding Equation 6 
(ај = [ај - [2 | 2.6] 8 


Now, substituting into Equation І, 
я I T | 
2EU = |P ia, | 6 FG (Pia; 
2 | | 30 
= [eja] н] [Р іа, сай 
El 15 а symmetric matrix. Its members are the coefficients of the 


strain energy terms when the energy is written as a function of applied 


and redundant loads. Partitioning [z] into four sub-matrices, 





k=; = 


„а | 22 гхғ 


Р = applied loads = 33 
r= redundant /oads=/6 


Rewriting Equation 29, 


2 AN 2 
глинени S/O ври а (31) 


Minimizing the strain energy with respect to the redundant, а’, yields 





ac TP) E 2PH, + 20,H,, (32) 
Rearranging { г, ў = — м] A G ў 
= EN [2 | => 


34 


The matrix of all the unknown internal loads is given by 


(а), = le) [Zee] tP) on 


TN 


баја, = 15] 12) 


S is the unit load distribution matrix for the indeterminate structure. 
Each Sin term is the value of the i'th member load when the. structure 
is subjected to a unit value of the m'th applied load. 

Using the notation of Reference 5, the strain energy can be 


written: 
ues | (35) 
И = 2 |А] c| (2] 
where C is the matrix of flexibility influence coefficients. Sub- 


stituting from Equation 11 into Equation 23, 


v =le 314 - 
therefore 
c= Z Е Gn 
= L |T EN] EFG ЕЯ 
= εν ШЇ 
and substituting EN ΠΝ EM from Equation 10, 
с = [М НА" | on 
Note that Equation 34 may be expanded to read 
а, = $44 Р, + 512 P, ο» Р (39) 
a, = 521 Р, T 522 Р, ша Son Es 


3⁄5 


Therefore, each internal load, a, , may be found by substituting a 
particular set of applied loads, {Р} . For any other loading condition, 
only the matrix of applied loads needs to be changed to obtain the new 
internal loads for that loading (flight) condition. This enables a 
rapid evaluation of the maximum load on any element over a range of 
applied loads. 

A computer program, originally set up by Reference 2 for the CDC 1604 
computer was rewritten in Fortran IV to fit the IBM 360/67 computer and 
appears in Appendix A. This program performs the complete matrix oper- 
ation outlined in Section C with input A and B from the flexibility matrix 
and the applied loads, P. Matrices А, В, F, and P appear in Tables IV, 


V, VI, and VTI. 
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EXPERIMENTAL RESULTS 


In order to ascertain whether or not the results obtained from the 
strain energy method were accurate for this type of thick skinned wing 
construction, it was necessary to have experimental data with which to 
compare the results of the idealized model. The wing is mounted in a 
jig which closely approximates the actual method of joining the wing to 
the aircraft fuselage, as shown in Figure 15. Since the wing is a 
continuous structure, the root fixity coefficient is not changed as when 
a single span wing is mounted on a completely rigid bulkhead. The 
loading condition investigated was a torque of 336,000 in. lbs. This 
load was considered to be well within the elastic range, but still of 
sufficient magnitude to give adequate strain levels throughout the 
structure for repeatable strain measurements. Some strain gages were 
installed on the wing when received from Ling-Temco-Vought Aerospace 
Corp., and others were installed at the Naval Postgraduate School. 
Figure 16 shows the wing with one cover skin removed for installation 
of strain gages. Wing section » = 34.5 was particularly well instru- 
mented, with strain gage rosettes placed on both sides of top and bottom 
skins to check for differential bending effects. Strain gage data were 
reduced to principal stresses and directions through computer program 
"Recrose", written in Fortran IV for the IBM 360/67 computer, as shown 
in Appendix B. The torque loads were applied with four hydraulic cylin- 
ders through a heavy fitting attached to the wing fold lugs. The applied 
torque was a e be perpendicular to the elastic axis. The load 
was applied at both wing fold axes; and to ensure a symmetrical load, 


the four hydraulic cylinders were supplied through a common manifold. 
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Load magnitude was determined by a manifold pressure gage in the 
hydraulic pressure system and checked by four Dillon dynamometers 
attached to each cylinder. The strain readings were taken with a 
Budd /Datran digital strain indicator (Model TC22). All gages were 
wired to a common junction panel as shown in Figure 17. All gages 
were balanced at a zero reading before applying the load. Eighteen 
gages were read at each loading in about one minute. This minimized 
the chance for drift in gage reading. All gages were rechecked for 
their zero reading after the load was released and sufficient time 
had elapsed for hysteresis effects to be minimized. All testing was 
done at night to eliminate the thermal effects of sunlight shining on 
the wing. Reference 2 contains all strain gage coordinates, gage fac- 
tors and resistances. It also shows that the strains are in the linear 
region by comparing the gage readings at a loading of 8000 psi (336,000 
in. lbs. of torque) with those taken at 7000 psi (294,000 in. lbs. of 
torque) and extrapolated to 8000 psi. 

Strain readings for the rosettes on the cover panels appear in 
Table VIII. The readings were averaged for the back-to-back gages 
and entered in program "Recrose" to compute the stresses. The normal 
and shearing stresses and the principal axes for the skin panels are 


also shown in Table VIII. 
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DISCUSSTON OF RESULTS 


There are two areas of comparison in this study. The first is the 
comparison of the theoretical results with the experimental results, and 
the second is the comparison of the theoretical analysis of the wing 
with and without including the center section. 

In the first area, the stresses as found in Table IX show fairly 
close agreement. In the panels where the best stress agreement occurs, 
such as panel 38, the strain readings on both sides of the skin are 
nearly equal. The strain readings of back-to-back gages were averaged 
and used to find the experimental stresses at each station. "Where the 
corresponding interior and exterior gages showed considerable уагіапсе 
in their strain readings, the experimental and theoretical stresses also 
had greater variance. The difference in strain readings is due to 
differential bending in the cover skin and possible torsional resistance 
of the pivot rib at the rear fitting, both of which were neglected in the 
theoretical analysis. The differential bending was neglected because 
the model used for the flexibility coefficient of the cover skin was a 
flat panel which carries only shear. The torsional resistance of the 
pivot rib was neglected because the effect of an unknown amount of slack 
movement about the wing hinge could not be incorporated into the 
theoretical analysis. The strain gages themselves are believed to be 
accurate to within 5 micro-inches of the true strain values. Strain 
readings found by Holgren and Comfort in 1963 were very close to those 
found in this study. The substitution of aluminum pads between the 
116 and the loading floor for the original plywood pads, as recommended 


in Reference 2, does not seem to have made much difference in the strain 
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values. The original gages were carefully installed and waterproofed 
which seems to account for their longevity. Most of the gages have 

been in place for over four years. The normal expected life-span of 
this type gage in industrial use is one or two years. The first factor 
mentioned, that of differential bending and torsional resistance of the 
pivot rib, seems most likely to be the cause of the difference between 
theoretical and experimental results. The differences do not seem to 

be too great since some of the references mention that errors of 100% 
are sometimes found in the case of thin skinned structures. More con- 
clusive comparisons could be made if the torque loading were increased 
and larger strain readings were available. However, it was decided that 
a significant increase in the torque loading beyond the 336,000 in. lbs. 
used in this study might permanently warp the wing and render it useless 
for further study. 

The second area of comparison is between the theoretical results 
found in this study, and those of Reference 2. Including the center 
section, from У, = 0 to Y. = 25.095 resulted in changing 25 flexibility 
influence coefficients and did not change the equilibrium equations. In 
effect, it moved the wing root from the pivot rib at Yes = 25.095 to the 
center rib at Y = 0. Most of the coefficients changed by 3 or 4%, and 
none of them changed by more than 137; yet the matrix of resultant loads, 
(Ат) ‚ Shows some considerable changes. In comparing (Ат) of Reference 
2 with the {Ат} matrix found in this study, the former loads will be 
designated as aj (Reference 2). The greatest change was in the load 
on the rear intermediate Беап, К.І.В. m (Reference 2) = 1010.4; and 


а,, = 10,263.5. This load is applied to element 24 and will influence 
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the stresses on panels 23 and 25. The flexibility coefficient associated 
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with this force has a much smaller change, 3 (Ref 72) 39555 


“аз” д 
апа 87 43 = 14.3406. The purpose of mentioning these changes is to 
show that the cumulative effect of minor changes in 25 out of 129 non- 
zero flexibility coefficients may be great. “This means that small errors 
in deciding upon the true area and length of the elements of the ideal- 
ized model may be significant. When the actual structure to be idealized 
has irregular cut-outs and joints, and reduction in skin, web, and: 

flange areas occurs in steps rather than línearly, considerable effort 
will be involved in just constructing an idealized model. By examining 
Table X, it is seen that inclusion of the center section relieved the 
stresses at the leading edge of the wing which was the area of greatest 
discrepancy in Reference 2. Both normal and shear stresses are much 
closer to the experimental values. The stresses in the rear panels are 
also improved, but it is felt that the pins at point B and C in Figure 2 
fit snugly enough to offer some torsional resistance of the pivot rib at 
the rear fitting, which was not accounted for in the theoretical analysis. 


Tables XI and XII list the resultant loads, a, through a, for the wing with 


9 
and without center section included. The agreement between theoretical 
and experimental stresses throughout the wing seems to bear out the val- 
idity of the matrix-force method,and the approach as described in this 
study in the prediction of stresses in a wing prior to its construction. 
For greater accuracy, it would be necessary to use more elements which 


would require the calculation of additional equilibrium equations and hence 


expansion of the matrices. 


CONCLUS IONS 


The matrix-force method of analysis using the strain energy theory, 
as proposed in Reference 1, is acceptable for a thick-skinned wing. 

It can be expected to predict the stresses in a complex, highly redun- 
dant structure, loaded in torsion, with a reasonable degree of accuracy. 
For increased accuracy, the long slender panels representing the cover 
skin should be broken up, and smaller elements used which would reduce 
the size of the structural grid but increase the number of unknowns. 
Reference 7 mentions a method of accomplishing this by dividing each 
panel into a number of narrow longitudinal strips which could be treated 
as rectangles with constant shear flow in each one. The difference in 
shear forces of adjacent strips would be taken by assuming concentrated 
stringer areas in the spanwise direction. This method would still 
neglect Poisson's effect. From the results of this study, it appears 
that a better panel would be one whose length and width were of similar 
size. 

In summation, it was found that the accuracy of the matrix-force 
method was improved considerably (compared with the results found in 
Reference 2) by the inclusion of 14 additional members which resulted 
in changes to 25 flexibility coefficients. Generally speaking, it must 
be possible to attain a very high degree of accuracy by fundamentally 
the same method as has been used in this study, by dividing the structure 
into a much larger number of elements. This involves a larger number of 
redundancies and much more numerical work. This could be programmed to 
be accomplished by the computer. An experienced stress analyst, who 


knows the stress distribution of most of the elements of the structure 
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on which he is working, would be able to concentrate the redundancies 
to the unknown points, thus reducing the computational work considerably. 
The experimental analysis would be improved by the addition of more strain 
gages on the inside of the wind cover skins. 

Because of the number of man - hours involved in removing and 
replacing the cover skins, however, it would not be practical to install 


additional gages. 
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THE EDGE SHEAR FORCES ON THE PANELS 


5.64 
6482 
7.44 
12606 
7.54 
7.06 
5.64 
6.52 
7.44 
7.66 
7.54 
7.06 
0.30 
14.3287 
13.5041 
1320895 
12.5754 
12% 2295 
11.7984 
11.8046 
12.1320 
L25456 
13.03% 
1225007 
14.3408 


TABLE III 


14.3408 
13205977 
#3 0807 
1245728 
1-2. νυ 
11.8046 
79.3021 
74.9206 
71.3846 
68.5890 
66.4567 
64.9038 


4.1535 
4.0130 
39220 
74.2206 | 3.7965 
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3.5868 





TABLE IV 
MATRIX A INPUT FORM 
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TABLE V 


MATRIX B INPUT FORM 
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TABLE VII 


MATRIX :P- 


Input ‘Form 


1500.00 
292,680 
714.00 
528.00 
109.00 

-210.50 

-2782.00 
-4523.56 
=5212.17 
-5568.74 
—5 32>. SG 
-4593.58 
-4654.94 
-4654.94 
-4593.58 
-5395.86 
-5568.74 
-5212.17 
-4523.56 

232 0.73 

4825.58 

53°5.02 

6 5.17 

6: 5.85 

3802.52 

DES. 21 

2320.73 

4825.58 

52395.02 

6125.17 

62 9.25 

5602.52 
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TABLE VIII 


EXPERIMENTAL RESULTS AT Y,, 
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. GAGE NO. ee AVERAGE — - 
| | O min 
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-922.41 
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26002.06 
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-576.74 
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7695.25 


457.66 


255937 


599.43 
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495.37 


= ЕТ, 62 
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TABLE IX 






Comparison of Results at Station Y,, = 34.5. Loading 







336,000 in. lbs. of torque. Principal axis is measured 








from load reference axis, 





the Center Intermediate Beam, 






which has 33° 58' sweep back. 





MAX. = MIN. MAXIMUM SHEAR 
NORMAL| STRESS STRESS 


Experi- | Theore- | Experi- Theore-| Experi-| Theore- 
mental tical mental tical mental tical 


15 1094.02 






PRINCIPAL AXIS 

























1605.47 | 1008.22 | 1333.88 


17 559009 A9 29087 
19 1047.07 895 57 
21 678.85 2 OMS 

23 -792.63 819.44 
25 070755 627.44 | 1033.12 














28 960.47 767.42 
38 457.66 626.57 512227 
32 813508 
34 599.43 551. 95 
36 =/ 67.00 
38 -1111.62 11052. 83 
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TABLE X 


COMPARISON OF THEORETICAL RESULTS AT STATION Y,, 


(1) WITHOUT INCLUDING CENTER SECTION, AS IN REF. 


(2) INCLUDING CENTER SECTION 


MAXIMUM MINIMUM .. MAXIMUM | 
PANEL | NORMAL STRESS -| NORMAL «STRESS ' - | `'SHEAR-STRESS | 
NO. 
la 









17 
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TABLE XI 
MATRIX OF ALL INTERNAL LOADS 


{Α1]-- [5] {Ρ} 
FROM IDEALIZED MODEL WITH CENTER SECTION INCLUDED 


174.54054260 


1325.45062773 26 266.9868164 1 
462. (32833008 27 1781.35253906 
474.22583308 28 3944 .20239254 
675.57104492 29 -84С7. 73046875 
1266.64038С86 30 -7190.99218750 
-597.06024414 31 -5754.24218750 
-3451.01147461 32 -4264.46093750 
~ 3498, 754150939 33 -2847.58456094 
-4256.8125C90C 3m 650.99536133 
- 5356. 39062 500 35 716.78784180 
-5181..39С6250 36 392.53369141 
-4892, 48437500 37 -136.54725647 
-7350.36718750 38 7% 1679.07592773 
-7291.8125690С 39 6767. 17968750 
-5829.66 796875 40 6785. 34375200 
-5258.51953125 41 5224.67187500 
-4711. 71093750 42 985.44262 695 
-4С13.54619141 43 10263.53515625 
-3421.25048828 44 3625.59106445 
197: 23242188 45 6134.81640625 
2 qo 57227 46 5193.66015625 
2 52125966231 дй 3034. 79296875 
-1755.8449797С6 . 48 4810.€0546875 
-5792.251227272 ^9 6030.46484375 
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TABLE XII 


MATRIX OF ALL INTERNAL LOADS 


{Ат} - [5] (РЈ 
FROM IDEALIZED MODEL WITH CENTER SECTION NOT INCLUDED 


1 -92,32380616 2έ = 11005 SO 1553 
2 1^92,21518555 27 -939.5744628< 
3 -60.23043£22 28 -859.61654336 
4 S707 C2 S62 198 26 -3524.745€0547 
5 -2t3.^8£3482C 30 -3361.27685547 
έ АЗС 31 -2227.84252<3< 
7 442 .2Е46ЗЕЕТ 32 -26517.405С263С 
е -2286.2754882€ 33 -2714.0425452] 
5 -2785.3С151367 34 767.75556С66 
re - 4119 .С97Е5Е2Е 35 504 .96675€€7 
11 4 163458432500 36 592.9165С361 
12 -4985.1992187* 37 -17.3076524€ 
13 -4569.667566765 38 -720.50$6222( 
14 ИАС, 27 228 ЕС 3¢ 8591.921875CC 
5 -4354,523431£0 4C 8726.14С62<СС 
16 -46Е86,596052765 41 7870.3671675С 
17 кы. оча 42 4969, 910312=С 
18 Bees. S 1 Cees 42 1010 .42R47€5€ 
С -4215.2639с626 44 5465. 73431 ЕСС 
2G -~ 2867. 74585252 45 8681.625СС0ССС 
ZH 5143.&6R750CCC 4€ а3257.22125СС< 
27 200.6325έ83έ 47 1980.82с3126С 
2 з -379.16625671 48 45C3.76171 E75 
24 А рет 4G 765.41543356 
25 -1254.2205£6€41 | 
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APPENDIX A 


PROGRAM ENERGY 


START 


CLEAR ARRAYS 
FORM А,В, ІС 
CALI GAUSS 3 


INVERT 4 , WRITE AT 


CAIL MATMPY 
CJ =A 


RETURN 


PRINT MATRIX 
SINGULAR 


В 





FORM G FRCM CJ ЗУ 
ADDING TRIVIAL RELATIONS 
WRITE G 


TRANSPOSE G 


CAIT MATMPY 
H= GF G, WRITE 


FORM Р 


PARTITION H 


ΕΙΝ 522 





“PITE H22 


CAII. GAUSS 3 


<-> YES 


"мү рт 1199 
1 - 


PRINT MATRIX 
SINGUI AR 
















, 18175; 


H221 











FORM H21 > 


CALL MATMPY 
ENP = H22I x H21 


= -ENP 







UTIN == I/EN 





CALL MATMPY 
S = G X UNITN 
WRITE S 







CALL MATMPY 
AR — EN XP 
WRITE AR 







CALL MATMPY 
AT — xP 
WRITE AT 







FORM C 


C=1/E ( H11+H12 XEN ) 
WRITE C 
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